Local Rigidity in Sandpile Models 
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We address the problem of the role of the concept of local rigidity in the family of sandpile 
systems. We define rigidity as the ratio between the critical energy and the amplitude of the 
external perturbation and we show, in the framework of the Dynamically Driven Renormalization 
Group (DDRG), that any finite value of the rigidity in a generalized sandpile model renormalizes 
to an infinite value at the fixed point, i.e. on a large scale. The fixed point value of the rigidity 
allows then for a non ambiguous distinction between sandpile-like systems and diffusive systems. 
Numerical simulations support our analytical results. 
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I. INTRODUCTION 

The concept of Self Organized Criticality (SOC)0Si 
has been proposed as a unifying theoretical framework to 
describe a vast class of driven systems that evolve "spon- 
taneously" to a stationary state characterized by power 
law distributions of dissipation events. While originally 
SOC systems were associated to the absence of tuning 
parameters, it is now cleamQ that most of the proposed 
systems show a critical behaviour only in the limit of an 
infinite separation of time scales, i.e. if some suitable 
parameter (e.g. -the dissipation and the driving rate in 
sandpile modelstj, the temperature in some growth mod- 
elsB, etc.) is set to zero. These results suggested a new 
definition of SOCEl as the theory of dynamical processes 
that lead a system to the critical steady state if the crit- 
ical value of the control parameter is zero. 

The stabilizing effect of a thresholds in the class of 
sandpile models is discussed by introducing a local rigid- 
ity defined as the ratio between the critical energy and 
the amplitude of the external perturbation, r = e c /Se. A 
finite rigidity allows the system to assume a large num- 
ber of metastable configurations and this let the .process 
show a power law distribution of avalanches. InEl it has 
been shown by means of numerical simulations that in 
the limit r — > the system becomes a diffusive one char- 
acterized by only infinite avalanches. In this paper we 
show how, in the framework of the so-calleii— Dynami- 
cally Driven Renormalization Group (DDRG)oEj, a finite 
rigidity in the microscopic dynamics is crucial in order 
for the system to reach spontaneously a critical station- 
ary state instead of a diffusive one. The basic idea is that 
at large scale the value of the rigidity does not depend 
on the microscopic value and that this "coarse grained" 
value allows for a non ambiguous distinction between a 
diffusive and a SOC system. 



II. THE RENORMALIZATION OF THE 
RIGIDITY 

In Sandpile Models we assign an integer variable ( "en- 
ergy") Si on each site j of a d-dimensional hypercubic 



lattice. At each time step an energy grain is added on 
a randomly chosen site, and then the system is allowed 
to relax according to a particular stability criterion de- 
pending on a fixed threshold (e.g. the energy, or the slope 
-defined as the difference between the energy of two near- 
est neighbours- exceeding a critical value). Infinite slow 
driving, i.e. an infinite separation of time scales, is built 
into the model: during the updating process the exter- 
nal input stops. The original BTW sandpile modeltl is 
a critical height model with e c — 4 and external input 
Se = 1. 

We consider a generalized sandpile model with criti- 
cal height e c and external perturbation Se for the micro- 
scopic dynamics. If e< > e c at time t then e, = at time 
t + 1 and the q nearest neighbours of site i increase their 
value of Ei/q. At the generic scale these quantities have 

the values ei^ and 8e^ k \ The height of the stable sites 
is less than Ec — Se^ k \ for the critical sites the energy 
lies between el — Se^ and el while for unstable sites 



it is larger than si ■ We study the behaviour of e y c' and 
fe( fc ) under the renormalization flow of the DDRG which 
could be formally written as: 



.(*) 



E (*+i) = /(4 fe ),fe( fe >;p«,pW) 



(1) 



where p^ (density of critical sites) and P( fc ) = 

iPi^ i f4^ > Ps^ i ) (Pn is the probability that in a top- 
pling the energy is distributed exactly to n nearest neigh- 
bours) are the parameters characterizing the state at the 
generic scale k. In principle we have to couple these two 
equations to those that determine the renormalization of 
p( fc ' and pW; for sake of simplicity, having extended the 
parameters space, we assume that the flow of the new 
parameters is orthogonal to the space defined by the old 
ones. We then accordingly fix and to their fixed 
point values p* and P* as computed in ref.[9]. Their val- 
ues are p* = 0.468 and P* = (0.240,0.442,0.261,0.057). 
The validity of this approximation will be checked by 
numerical simulations. Our transformation will have the 
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FIG. 1: Critical, stable and unstable sites. 
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FIG. 2: The a configurations with relative statistical weights. 



form: 



fe( fc+1 > = .g(4 fc ),fcW;p*,P*) 



(2) 



The symbols used are shown in Fig. [l] and |2|. 

Let us consider now a single path as an example of 
the transformation we want to obtain. We have a 2 x 2 
cell characterized by a label a specifying the number of 
critical sites. The a — 1 configuration is not considered 
because of a spanning rule that imposes the exclusion of 
those processes that do not relax within the cell before 
they transfer energy outside. Every configuration a has 
a mean height (e a ) that we calculate in the next section. 
One of the possible evolution from the configuration a = 
3a state is shown in Fig. |^. 

This means that at the scale k + 1 the mean height of 
the considered cell is critical with probability 



— (a) 
7T = tlT - 



(3) 



As a matter of fact, the cell has a probability ir to be a 
critical site at the scale k + 1. The complete renormal- 
ization flow is therefore a sum over all the configuration 

and over all the possible paths of terms like ir ■ (e^l 3 ). In 
the same process we observe that with probability ir the 
quantity given outside of the cell, fe( fc+1 ', is 3/2 • fe( fc ). 
3 flows of energy 8e^ are lost outside the cell passing 
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FIG. 3: One of the possible time evolutions starting from an 
a — 3 configuration and transferring energy outside of the 
cell; for every step the relative probability is indicated. 



through 2 sides. It means that for this transformation 
we have to distinguish paths transferring the same en- 
ergy through a different number of sides. The general 
transformations are 



r (fc+i) 



(4) 



£ — ' m 



where Pi iCC are the probabilities relative to time evolu- 
tions that transfer energy outside of the cell, while P™^ m 
refers to a process that let n quanta going through m 
sides. 

We compute now the average energy per site in all 
the configurations. We generalize the simple mean field 
argument described in ref.[ll], making the simpler non 
trivial assumption for the probability distribution of the 
energy: P(e) — a + be. The normalization 



deP{e) 



1 



(5) 



implies that b = 2(1 — ae c )/e^. 

The average energy is then expressed in term of the 
parameter a: 

{e)=J deeP(e) = ^e c -^ae 2 c . (6) 

Following ref.[ll] we find that a ~ 0.277/e c . The average 
energy of stable and critical sites is then: 
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FIG. 4: Time t: the probability relative to the paths in figure 

is |pi + 5P2 + jP3 +P4- 
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If we put a = 0.277/e c in these expressions we obtain 

(e) stable = 0.6205 Ee - 0.539 <5e (9) 



(8) 



FIG. 5: Time t + 1: the probability associated with the paths 
in figure is §pi + |p2 + f>3 + P4- 
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(10) 



The average energy per site for different configurations 
can now be calculated as 

(e a = 2 ) = | [2 (e) cr mcal + 2 (e) stable] = 0.810 e c - 0.359 Se 

(11) 



(£ Q=3 > = J [3 {e)critical + (s) stable] = 0.905 £ c - 0.269 fc 

(12) 

(e a=4 ) = (e) C riticai =s c - 0.1785 Se (13) 
Equations for the renormalization flow have the form 



8e {h+x) = c5e (k) 



(14) 



where the coefficients a, b and c obviously depend on the 
fixed point parameters. If we divide the two equations, 
recalling the definition of the rigidity, we obtain 



r (fc+i) = ° r (k) _ i 
c c 



(15) 



The fixed point of the rigidity corresponds to setting 
r (fc+i) = r {k) _ r *. 



(16) 



a — c 



Next section deals with the calculation of the coefficients 
a, b and c. 

III. ANALYTICAL RESULTS 

To calculate the coefficients a, b and c we have to distin- 
guish between 4 different kind of configurations, each one 
with the relative weights w a (p) as calculated in ref. [9] . 

• a = 2 



We consider at time t one of the toppling represented 
in Fig. H and possible developments at time t+ 1 in Fig. ||. 

All the paths described transfer energy outside the cell, 
the associated probabilities are given in captions. Con- 
figuration a = 2 then contributes with the term 



w a=2 (p) ■ (ei1 2 ) • ( 7P1 + t;P2 + ~P3 +P4 ) x 



1 



(17) 



X ' 2 Pl 6 P2 +P3+Pi 



to the renormalization of the critical energy. Moreover, 
each of the described paths transfer only a quantum of 
energy outside of the cell, so 

fe(fe+i) f\ 1 3 \ 

= w a=2 (p) ■ f ~ Pl + -p 2 + -p 3+P4 \x 



(18) 



X ' 2 Pl 6 P2 +P3+Pi 



The details of the calculations regarding paths from 
configurations a = 3a, 3b and 4 are given in the ap- 
pendix. 

The final result is obtained by summing all these con- 
tributes that we have explicitly calculated with the re- 
spective weights. We find for the coefficients the follow- 
ing values: 



= 0.26 ±0.03 b = 0.10 ±0.02 



0.26 ±0.05 (19) 



We finally obtain that a and c are equal, even though 
the errors associated are relatively large. These errors 
are related to the parameters of the fixed point. Since 
they have only 3 significant digits, that results an error 
of order 10~ 3 . 

We then find that the fixed point spontaneously 
reached by a sandpile model corresponds to an infinite 
value of the rigidity. It means that every finite value for 
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FIG. 6: The behaviour of the critical energy at the scale 
k in the BTW sandpile model. 



FIG. 7: Numerical renormalization of Se in the BTW. 



the microscopic dynamics renormalizes to an infinite one 
on a large scale. 

The difference with a diffusive system is now clear: in 
the diffusive case a particle added from outside can al- 
ways goes out of the system from the boundaries, while in 
a SOC system this is avoided by the value of the rigidity 
as bigger as observed from a large scale. 
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IV. NUMERICAL RESULTS 

We_present hece-the results of simulations made on 
BTWld and Zhango models. We have measured for these 
systems the critical energy and the average flow trans- 
ferred from the boundaries at a generic scale k. Critical 
energy is defined as follows: we consider an external per- 
turbation on a site in a coarse grained cell (fc-cell) of size 
iWxlW and we look at the avalanches starting from this 
coarse-grained cell. If the avalanche transfers energy out- 
side the given fc-cell we compute the average energy Be 
of that fc-cell before the toppling; as a matter of fact, that 
fc-cell is critical with a probability given by the frequency 
of the avalanches transferring energy outside the fc-cell it- 
self. Moreover, we define the energy transferred at 
the scale as the number of grains going outside the 
fc-cell divided by the number of boundaries crossed and 
by the length (this is the trivial scaling). 

If we look at numerical results (Fig. ||) we note that 
the critical energy at large scale is the average energycHl 
of the system and this is not surprising since at the scale 
of the system size this is exactly what we compute. As 
for the behaviour of Se (Fig. |?]) one gets that in the limit 
k — » oo it goes to zero. This result confirms the idea that 
the fixed point value for the rigidity is infinite, signalling 
a clear difference with respect to usual diffusive systems. 
The same qualitative results are obtained for the Zhang 
modeO (See Fig. | and |). 
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FIG. 8: The behaviour of the critical energy at a generic scale 
k for the Zhang model. 



V. CONCLUSIONS 



In summary we have investigated, both in the frame- 
work of the Dynamically Driven Renormalizaton Group 
scheme and by numerical simulations, the role of the so- 
called local rigidity in sandpile models. The local rigid- 
ity is defined as the ratio between the critical energy and 
the amplitude of the external perturbation, both at the 
microscopic scale. It turns out that, under the renor- 
malization flow equations of the DDRG, the local rigid- 
ity renormalizes to an infinite fixed point value both for 
the BTW model and the Zhang model, while for a typi- 
cal diffusive system one would expect a vanishing value. 
Numerical simulations of coarse-grained sandpile models 
confirm this picture. 

Authors acknowledge EU contract FMRXCT980183 
and G.C. acknowledge FET Open project COSIN IST- 
2001-33555. 
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FIG. 9: Numerical renormalization of Se in the Zhang model. 
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VI. APPENDIX 



• a — 3a 

In Fig. [H^ we show all the possible paths starting from 
the perturbation of the critical site. 

The probability relative to the first event at time t is 
the same as before, ^pi + \p 2 + jP3 + P4- At time t+1 
we have a) \p\ + ^p 2 for the first two events in figure, 
\p2 + jP3 for the third (they are distinguishable only 
with respect to renormalization of Se); b) ip2 + jP3 for 



the first two events, jp^ + P4 for the third; 



\P2- 
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At time t + 2 we have \p\ + \p 2 + \p-$ for the first two 
and + \p3 +P4 for the third; in the renormalization of 
e c we have to consider only the sum ^p x + |p 2 + P3 + P4- 
So we can write 



,(k+i) 



• J rW ( S S (p)-{e^l s )-( - Pl + ^p 2 + |p 3 +P4 } 



.(*) 



\pi + \P2 +P3 +P4 J U + (\pi + \pi 
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FIG. 10: Time t: the probability relative to the paths in figure 
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P2 + P3 + P4 



(21) 



for the contributions coming from the a = 3a configura- 
tion. 

• a = 3b 

We can compute the probabilities of the shown in 
Fig. 0. 

At time £ we have to distinguish between two kinds 
of events: the toppling on a single site (this term is 
considered two times to include the toppling on the 
other site) and the one on two sites. The statistical 
weights are (\p x + \p 2 + jPs) and (\p 2 + \pz + Pa) re- 
spectively. At time t + 1 both the events leave one flow 
going outside the cell, so we have: 
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FIG. 11: Evolutions for the a = 36 configuration. 
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FIG. 12: Events a and b for the a — 4 configuration. 
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(23) 
• a = 4 

In a cell with all critical sites the number of possible 
time evolutions grows enormously. In Fig. [l2| we show 
the first two possible events, with probabilities p{a) = 
\P\ + \P2 + jP3 and p(b) = ±p 2 + hp3 + Pa- The paths 
coming from a are shown in Fig. O and the relatives 
probabilities in Tabic |[ 

Now we have to compute 224 paths from the b event. 
The final result, avoiding details, for the renormalization 
of the critical height is 
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FIG. 13: Possible paths from event a. 
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TABLE I: Probabilities relative to the events in Fig. 
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If we introduce the quantities 
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then we can write down the renormalization of Se as 

§E(k) = ■■■+ W a =i{p) ■ {2 7T 5 [7T 3 + 
7T4 ^2 7T 7 + 7T 4 ^1 + 7T 2 + ^7T 5 + ^7T 6 ^ ^ + 

+7T"9 ^7T 5 + |7T 6 ^ + 7T 7 7T 2 + 7T 7 (2 + 1) + 

7T, | - + \ \ ■ \ \ ■■■■■ I - «i !! - - ; ! - I + 
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(25) 
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